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Abstract We present timing measurements, astrometry, and 
high-resolution spectra of a number of nearby, thermally emit- 
ting, isolated neutron stars. We use these to infer magnetic 
field strengths and distances, but also encounter a number 
of puzzles. We discuss three specific ones in detail: (i) For 
RX J0720.4-3125 and RX Jl 308.6+2127, the characteris- 
tic ages are in excess of 1 Myr, while their temperatures 
and kinematic ages indicate that they are much younger; 
(ii) For RX J1856.5— 3754, the brightness temperature for 
the optical emission is in excess of that measured at X-ray 
wavelengths for reasonable neutron-star radii; (iii) For RX 
J0720.4— 3125, the spectrum changed from an initially fea- 
tureless state to one with an absorption feature, yet there was 
only a relatively small change in T e g. Furthermore, we at- 
tempt to see whether the spectra of all seven sourced, in six 
of which absorption features have now been found, can be 
understood in the context of strongly magnetised hydrogen 
atmospheres. We find that the energies of the absorption fea- 
tures can be reproduced, but that the featureless spectra of 
some sources, especially the Wien-like high-energy tails, re- 
main puzzling. 
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1 Introduction 



One o f the great benefits of the ROSAT All-Sky Survey dVoges etaO 
1996) is that is has provides an unbiased sample of all classes 
of nearby neutron stars (limited only by their age and distri- 
bution of the local interstellar medium). Particularly inter- 
esting is the discovery of the group of seven nearby, ther- 
mally emitting, isolated neutron stars (INS; for a review, see 
Haberl, these proceedings). 

The INS form the majority among the ne arby neutron 
stars (t ypical distances are less than ^500 pc; Kapl an et alJ 
2002b; see also Posselt, Popov, these proceedings), yet are 
atypical of the neutron-star population represented by radio 
surveys: while pulsars detected by their thermal emission 
all have normal periods of less than a second, five out of 
the seven INS have periods about ten times longer (the re- 
maining two appear to have no pulsations despite inten sive 
searches: iRansom et alJl2002t Ivan Kerkwii k et alJl2004h . A 
number of models — accretors ilWanglll997h. middle-age d 
magnetars dHevl & Kulkarnill998||Hevl & Hernauistl l999). 
long-period pulsars ( Ka plan et al J 12002a : IZane et alll2002l) 
— have been suggested to explain these objects. 

A prime reason for studying the INS is the hope of con- 
straining fundamental physics at very high densities: neu- 
tron sta rs are natur al laboratories for quantum chromody- 
namics (Rho 2000). The overall goal is to determine the 
masses and radii of a number of neutron stars and hence 
con strain the equation of sta te (EOS) of ultra-dense mat- 
ter dLattimer & Prakashl20 00: Lattimer, these proceedings). 
For the majority of known neutron stars (i.e., radio pulsars), 
this is complicated by the non-thermal emission that domi- 
nates the spectrum, but for the INS this is not the case: the 
X-ray spectra show thermal emission only. Hence, much ef- 
fort has been spent trying to derive constraints from the INS 
(iRurwitz et ahkoOlLl^lPrlke et al fcOOfllPons et al J2002h. 
The constraints have not been very meaningful, however, be- 
cause the data could not be interpreted properl y: they just 
do not fit any c urrent realistic models (M otch et al]l2003t 
ane etall2f 
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To make progress in understanding the thermal emission, 
we need first to know the basic ingredients: the elemen- 
tal abundances, the temperature distribution, and the mag- 
netic field strength. Furthermore, to use the thermal emission 
to infer radii, we need information about the distance. For- 
tunately, observational clues are now becoming available: 
broad absorption features at energies of 0.3-0.7 ke V have 
been discover e d in the spectra of six of the seven INS ( Haberl et 
120031 l2004bfl Ivan Kerkwijk et alJl200l IZane et alJl2005t 
see Haberl, these proceedings), and, as described below, mag- 
netic field strengths have been inferred from timing solutions 
and new or improved parallaxes have been measured. 

The outline of this contribution is as follows. First, in 
|J2| we present timing solutions for RX J0720.4— 3125 and 
RX J1308.6+2127, and discuss the resulting estimates of 
the magnetic field strengths and characteristic ages. Next, in 
fj3l we describe new parallax distance measurements for RX 
J1856.5-3754 and RX J0720.4-3125. For the former, these 
resolve previous conflicting results, but also raise a puzzle: 
a rather large radius or high brightness temperature inferred 
for the optical emission. In i|4] we turn to high-resolution X- 
ray spectra, comparing spectra of RX J0720.4— 3 125, before 
and after its spectral change, with those of RX J1308. 6+2127. 
In we attempt to interpret the observations assuming the 
sources have gaseous atmospheres, focussing on hydrogen, 
but also briefly discussing the possibility of helium. We sum- 
marise and discuss future work in Sj6] 

From here on, we will refer to the INS in the text using 
abbreviated names: J0420 for RX J0420.0-5022, J0720 for 
RX J0720.4-3125, J0806 for RX J0806.4-4123, J1308 for 
RXJ1308.6+2127 =RBS 1223, J1605 for RX J 1605. 3+3249, 
J1856 for RX J1856.5-3754, J2143 for RX J2143.0+0654 
= RBS 1774 = RXS J214303.7+065419. 




Period (s) 



Fig. 1 P-P diagram, showing radio pulsars (points) and magnetars (di- 
amonds); selected objects are labeled. Also shown are the five INS with 
periodicities: RX J0720.4-3125 and RX J1308.6+2127 are shown 
by the stars, while RX J0420.0-5022, RX J0806.4-4123, and RX 
J2143. 0+0654 are the arrows at the top (since P is unknown). The di- 
agonal lines show loci of constant dipole magnetic field and spin-down 
age, as labeled. 



2 Timing solutions 

Until recently, the typical magnetic field strength of the INS 
was just a guess, with a wide range of possibilities (10 10 - 
10 15 G) - similar to the wide uncertainty in magnetic field 
strength implied for the different demographic models of the 
INS, although based on the long spin periods fields of a few 
10 13 G were considered most likely. 

We have improved upon this situation using dedicated 
timing observations with Chandra, which, combined with 
Chandra, XMM-Newton, and ROSAT archival observations, 
allowed us to determine phase-coherent timing solu tions for 
J0720 and J1308 stretching back at least 5 years ( Ka plan & van 
2005 albl) . Using the periods and the period derivatives to 
place these objects on the classic P-P diagram (Fig.0, one 
sees that they are intermediate between radio pulsars and 
magnetars (in line with the idea that they are long-period 
pulsa rs whose radio beams do not c ross our line of sight; 
Kapla n et all2002atlZane et all2002l) . 

From the solutions, assuming the sources spin down by 
magnetic dipole radiation, one infers that they have simi- 
lar magnetic fields: B = 2.4 x 10 13 G and 3.4 x 10 13 G for 



J0720 and J 1308, respectively. As will become clear below 
(a3, this agrees quite well with the magnetic field strengths 
inferred from the absorption lines. 

One also infers characteristic ages: T c = P/2P = 1 .9 Myr 
for J0720 and 1.5 Myr for J1308. These are puzzling, since 
they are substantially longer than expected based on cool- 
ing: From standard cooling curves ( Page et al. 2004[), the 
observed temperatures around 90 eV (10 6 K) correspond to 
ages of a few 10 5 yr. Even if one takes into account that 
the black-body temperature likely overestimates the effec- 
tive temperature - as is clear from the fact that the extrapo- 
lation of a black-body fit to the X-ray data underpredicts the 
optical - one is very hard-pressed to find an age in excess of 
10 6 yr; at 1.5 Myr, the effective temperature should be below 
20 eV. 

For J0720, there is an a dditional age estimate from kine- 
matics dMotch et alJl2003b lKaplanll2004l also Motch, these 
proceedings): tracing back its proper motion, the most natu- 
ral birthplace is the Trumpler 10 association; it would have 
left about 7 x 10 5 yr ago. (Note that a possible origin in the 
Scorpius O B associations about 1 .5 Myr ago has also been 
Kmfeslft , but the new parallax and proper motion we de- 
rived [partly described in below] make this less likely. 
Furthermore, for J 1 856, which is cooler than J0720, the kine- 
matic age of about 4 x 10 5 yr is not in doubt.) 

Of course, the above discrepancy may simply mean that 
the characteristic age is a poor estimate of the true age. For 
breaking with V <x v", where V is the spin frequency and n 
the so-called braking index (equal to 3 for magnetic dipole 
radiation), the true age is f = (1 - P Q / P) n ~ l (P / (n - \)P). 
Thus, one can obtain ages («T either if the initial spin pe- 
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Fig. 2 Phase residuals for timing measurements for RX J0720.4— 3125 
relati ve to the best-fit quadrati c (v = ) model. The data include those 
from Kaplan & van Kerkwiikf <2005aft . plus a number of additional 
points from Chandra with LETG/HRC and one additional point from 
XMM-Newton with EPIC-PN (all at MJD > 53500); other EPIC-PN 
data exist (Haberl, these proceedings) but they are not yet public. 
The fit has %v = 11.93 and the the residuals have rms = 0.39 s. We 
also show alternate fits: the best-fit cubic model (v ^ 0; dotted, black 
curve), which has rms = 0.36 s and xl = 5.12; the best-fit periodic 
model (dot-dashed, blue), which has a period of 4.3 yr and rms = 0.23 s 
and x 2 = 2.35; and a simple glitch model (dashed, red) with the glitch 
occuring at MJD 52821, which does not substantially improve the fit 
over the quadratic model. Note that there appears to be high-frequency 
residuals that are not fit by any of these models, as can be seen by the 
disagreement between simultaneous XMM observations with different 
instruments. Whether these are due to the energy dependence of the 
pulse profile or some instrumental effect, we cannot say, but it sug- 
gests that caution should be taken in interpreting the timing residuals. 



riod Pq is close to the current one (i.e., the neutron star was 
born spinning slowly, but in no other systems is there evi- 
dence for Po > 1 s), or if n is substantially larger than 3. Lyne 
(these proceedings) presented evidence for values of n sub- 
stantially different from 3, although most values were less, 
implying characteristic ages that are longer than the true age, 
contrary to what is required here. 

As an alternative, we noted that one way of obtaining 
Po — P, would be to have the neutron star undergo a phase 
in which it was accreting, either from a companion (which 
later disappeared, e.g., in a supernova explosion) or perhaps 
a fall-back di sk such as that di scovered around the AXP 
4U 0142+61 (W ang et alJl2006h . Intriguingly, the equilib- 
rium spin period, P eq « 5s(B/10 13 G) 6 / 7 (M/M Edd )" 3 / 7 , is 
roughly equal to the current observed periods for a magnetic 
field of a few 10 13 G and an accretion rate M close to the 
Eddington rate MecM- 

Finally, comparing the timing residuals, we find that for 
J0720, they are ^-0.3 s, far larger than the measurement er- 
rors, while for J1308, they are consistent with the measure- 
ment errors, at ^-0.01 s. The larger residuals for J0720 have 
been ascribed to precession ( Habe rl et alJ2006t also Haberl, 
these proceedings). To verify this, we tried including differ- 
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Fig. 3 Improvement in Xv f° r a timing model for RX J0720.4— 3125 
that includes a periodic component, as a function of period, fit to the 
data shown in Fig. [2] The Xv (solid, black curve) has its deepest mini- 
mum around 1580 days (4.3 yr); there may be a second minimum near 
3000 days (8.2 yr), close to the ~7 yr period suggested by Haberl et 
al. (2006), but with these data it is clearly less significant. Also shown 
is the fractional change in v for the different fits (dashed, red curve); 
it changes by < 2% for the full range of periods considered here, and 
thus does not influence the conclusions regarding the magnetic field or 
characteristic age. 



ent terms in our timing model (see Fig. and we indeed 
find that adding a periodic component improves the fit dras- 
tically, with the reduced % 2 decreasing from %l = 11.9 to 
2.4. Trying different periods, however, the best period ap- 
pears to be 4.3 yr (Fig. [5}, and not ~7 yr as inferred from the 
spectral changes (at 7 yr, the fit is better than quadratic, but 
not that much different from a higher-order polynomial). We 
caution, though, that the timing residuals shown by J0720 
are not exceptional: they are in line with trends seen for radio 
pulsars and similar apparent periodicities can be seen in the 
residuals of some of the Anomalous X-ray Pulsars (Kaspi, 
these proceedings). 



3 Parallax Measurements 

A parallax measure ment for the brightest INS, J1856, was 
first attempted by IWalterl d200 lh . using three observations 
with the Planetary Camera onboard the Hubble Space Tele- 
scope (HST); the resulting parallax implied a distance of 
^60 pc. The measurement was tricky, however, and a much 
larger dist ance of 140 pc was derived from the same obser- 
vations by Kap lan et al.ld2 002b): a lar ger dis tance, of 1 17 pc, 
was also found by Walter & Lattimen (12002 ). who redid their 
analysis and included a fourth PC observation. 

In order to obtain more accurate distances, we have used 
the High-Resolution Camera (HRC) of the Advanced Cam- 
era for Surveys (ACS) onboard HST. This camera is more 
sensitive than the PC and has a smaller pixel scale, so that 
undersampling of the point-spread function is much less of 
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an issue. Furthermore, across the pixels, the sensitivity is 
more uniform, reducing the variability in the point-spread 
function with pixel phase; as a result, m uch more accurate 
astrom etry can be done with ACS/HRC (lAnderson & Kina 
2006). We obtained images of J1856 and J0720 in the blue 
F475W band, visiting each source eight times over two years. 

For J 1856, our analysis of the HRC data is virtually com- 
plete (Kaplan, van Kerkwijk, & Anderson, in preparation). 
In order to obtain as accurate a parallax as possible, we have 
taken into account the parallactic motion of the background 
reference stars, by determining their photometric parallaxes 
(assuming they are main-sequence stars; for the less dis- 
tant stars - generally the brighter ones with strong weight 
- we confirm the photometric parallaxes astrometrically). 
With that, from the HRC data alone, we determine a par- 
allax n = 6.2 ± 0.6 mas, corresponding to a distance d = 
161 pc. We are currently trying to improve the measure- 
ment further by including the PC data. 

For J0720, a first analysis of the HRC data has just been 
completed. For this source, parallaxes of the background 
stars are much less important, since it is at low Galactic 
latitude and most objects are distant. Our preliminary par- 
allax is n — 3.0 ± 1 .0 mas, corresponding to a distance d — 
330±™ pc. 

The factor two ratio in the distance s to J 1856 and J0720 
is consistent with what was expected by Kap lan et alJ (12002b) 
under the zeroth-order assumption that the optical flux for 
different sources scales as °= T(R/d) 2 , that the radii R are 
similar, and that the temperature T in the region of the atmo- 
sphere emitting the optical emission scales with the temper- 
ature determined from fits to the X-ray spectrum. The dis- 
tances also compare well with the distances of 135 ±25 and 
255 ± 25 pc inferred from the run of H I column density with 
distance (Posselt, these proceedings). 

With our distances, we can estimate the radii for the 
two sources. We start by simply using the black-body fit 
to the X-ray spectra. For J 1856, one finds kT = 63 eV and 
Roo/d = 0.0364 kmpc -1 , which, with our new distance, im- 
plies a radiation radius Roo ~ 6.5 km. This is smaller than a 
typical radius of a neutron star, but this is not unexpected, for 
two reasons. First, for the most likely atmospheric compo- 
sitions, the opacity decreases with increasing frequency. As 
a result, at X-ray energies one sees relatively hot layers and 
a fit to the X-ray spectrum will thus overestimate the effec- 
tive te mperature and underestimate the radius dPavlov et alJ 
1996). Second, the temperature distribution likely is not uni- 
form, in which case the area inferred from the X-ray emis- 
sion would simply correspond to that of the hotter parts. 

In the above picture, one expects the optical emission 
to be in excess of the extrapolation from the black-body fit, 
since it arises from a cooler layer and from a larger area. 
And indeed, the spectral energy distribution, shown in Fig.0] 
shows an excess. It poses a possible problem, however, since 
the optical excess is a factor 7, which implies a radiation 
radius of 7?oo. pt = 17(7x/r opt ) 2 km (where T opt and Tx are 
suitable averages of the temperatures of the optical and X- 
ray emitting regions, respectively). Given that one expects 



?opt < Tx, the optical emission thus seems to imply that the 
radiation radius is quite a bit larger than 17 km. Yet, for most 
reasonable equations of state, a typical neutron star will ha ve 
a smaller radiation radius (e.g.. lLattimer & Prakashl2 001). 

Of course, the above discrepancy may simply reflect our 
lack of understanding of neutron star atmospheres in strong 
magnetic fields: the temperature of the optical emission re- 
gion may be larger than expected. In order to see what would 
be required, one can reverse the process: assume that the 
neutron star has a 'standard' mass and radius, and calculate 
the brightness temperature at each energy assuming that the 
emission originates from the whole surface. In Fig. [5] we 
show the result for Roo = 14.7 km (which is the value one 
obtains for M = 1 .35 M e and R = 12 km). We see that this 
confirms the above reasoning: in order to produce the opti- 
cal emission, the temperature in the emission region has to 
exceed 70 eV, i.e., be higher than that in the X-ray emitting 
region. 

For J0720, a fit to the X-ray spectrum from 2000 (i.e., be- 
fore the appearance of an absorption line) gives kT = 85.7 eV 
and Roo/d = 0.0170 km pc" 1 . Taking the distance at face 
value, the implied radiation radius is 5.7 km, a little smaller 
than that of J 1856, but easily consistent within the 30% un- 
certainty due to the parallax measurement error. Since the 
optical excess is similar, the radiation radius for the opti- 
cal emission is again large. In this case, howev er, the o ptical 
emission does not follow a Rayleigh-Jeans tail ( Kapla netal] 
120031 : iMotch et all2003l: Fig.||}, and hence it is not clear that 
the emission is from the surface. This can also be seen from 
the brightness temperatures (Fig.|5j, which is not constant in 
the optic al/ultraviolet range. 

A further puzzle raised in comparing the sources, is that 
despite the fact that the X-ray emission areas are rather sim- 
ilar, the pulsation properties are very different: J0720 shows 
clear pulsations, with a pulsed fraction of 1 1% (lHaberl et alJ 
while J 1856 shows no pulsations, to a limit of ~ 1% 
(Rans om et alJ l2002: Bu rwitz et alJl2003h . This may reflect 
differences in geometry; for isotropic emission from two op- 
posite magnetic poles, there is a fair range in parameters for 
which no pulsations would be observed (e.g., Belob orodovl 
2002). Of course, the presence of the pulsations constitutes a 
warning about the brightness temperatures shown in Fig. |5] 
if the X-ray emission does not arise from the whole surface, 
the true temperatures will be higher than those shown. 



4 LETG spectra 

The study of the X-ray spectra of the INS has made great 
strides with the advent of Chandra and XMM-Newton. Both 
the CCD instruments, in particular EPIC-PN, and the grating 
spectrometers LETG and RGS have been used extensively. 
Here, we focus on the grating instruments (for the exciting 
results from EPIC-PN, see Haberl, these proceedings). We 
will only discuss results from LETG, since that instrument 
covers the full range of energies at which INS emit and since 
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Fig. 4 Spectral energy distributions for RX J1856. 5-3754 (left) and RX J0720.4-3125 (right; with the X-ray spectrum from before the ap- 
pearance of an absorption feature). For both, the X-ray points are from LETG spectra, the dark blue points from HST, and the cyan points from 
ground-based observations. The optical and ultraviolet spectra for RX J1856.5— 3754 are from VLT and HST, respectively. The black, drawn 
curves represent the best-fit black-body models to the X-ray data; the dotted curves are the same model without interstellar extinction. 





E (eV) 



E (eV) 



Fig. 5 Brightness temperatures for RX J1856.5— 3754 (left) and RX J0720.4— 3125 (right), assuming our parallax measurements are correct and 
that the emission arises from a neutron star with radiation radius R„, = 14.7 km (which is the value one obtains for M = 1 .35 M,-.-, and R = 12 km). 
The symbols and colours are as in Fig. |4] (the model curves are not at a constant temperature since the best-fit radiation radius is not equal to 
14.7 km). One sees that the optical emission requires temperatures at least equal to those required for the X-ray emission. Note, however, that 
for RX J0720.4— 3125 it is not clear the emission is thermal. 



the calibration of the RGS at longer wavelengths has been 
rather problematic. 1 

So far, three sources have been observed with LETG. 
By far the best spectrum is of J 1856, taken using 500 ks of 

1 The problems appear to be largely solved with the 2006 June 30 
release of the Scientific Analysis Software. 



director's discretionary time. Unfortunately, and puzzlingly, 
the spectrum appears completely featureless, a nd is well de- 
scribed by a black-body model dBurwitz et alJl200lL l2003t : 
iDrake et alJ2002l : lBraie & Romanil2002l) . The second bright- 
est source, J0720 has been studi ed extensively as w ell. A 
first spectrum was taken in 2000 iKaplan et 

1012003), 

when 
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(A) 



Fig. 6 Spectra of RX J0720.4-3125 and J1308.6+2127 taken with 
LETG. The magenta points connected with the line are the average 
of spectra of RX J0720.4-3 125 taken after 2004 (about 300 ks). Com- 
pared to the earlier spectrum (red points), the spectrum is harder and 
has developed an absorption feature. On the other hand, compared to 
RX J1308. 6+2127 (blue points, multiplied by two in order to match at 
the short-wavelength side), the feature is rather weak, even though the 
inferred temperature is similar (as is also clear from the good match at 
short wavelengths), and so is the magnetic field inferred from timing. 



the sp ectrum was featureless, and a second in 2004 (IVink et al.l 
2004]), t o confirm the change in spectrum discovered with 
RGS bv lde Vries et al.l (12004). During 2005, the source was 
regularly monitored by us, for 300 ks total, in order to look 
for further spectral changes and to see if more than one spec- 
tral feature was present. Finally, a lOOks observation was 
made of J1308 in guaranteed time by the MPA group. 

In Fig.|6] we show the LETG spectra of J0720 and J1308. 
For the former, we show separately the 'before' spectrum 
(from 2000, before the appearance of an absorption feature) 
and the 'after' spectrum (the average of all spectra taken af- 
ter the change; within our statistics, the individual spectra do 
not differ). 

The comparison of the three spectra raises a number of 
questions. First, for J0720, how can a relatively small change 
in temperature (from ~ 80 to 90 eV) cause the appearance 
of a pronounced absorption feature? Second, independent 
of the state of J0720, why is the absorption feature much 
less strong than that in J1308, despite the fact that their tem- 
peratures and magnetic field strengths (from timing, are 
rather similar? 

Considering first the change in J0720, the simplest possi- 
bility would be that the increase in temperature corresponds 
to a relatively large change in the ionisation and/or dissocia- 
tion balance. If so, this might give a clue to the correspond- 
ing energies of the matter in the atmosphere. The alternative 
would be that the region that heated up (or appeared in view) 



has either a different composition or a greatly different mag- 
netic field strength compared to the regions that dominated 
the spectrum before the change. Neither possibility seems 
particularly appealing. 

Whatever the physical reason for the appearance of the 
absorption feature, another more basic question is whether 
the change corresponds to a global change in the properties, 
or whether, instead, only a fraction of the surface changed or 
if ones viewpoint changed. A global change might occur if, 
e.g., heat was deposited deep inside the neutron star. In con- 
trast, a change in a limited area would be expected if heat 
were deposited near the surface (with the affected area per- 
haps increasing in size with time), or if the so urce were pro- 
cessing and different reg ions came into view dde Vries et alJ 
2004; Haberl et al. 2006; Haberl, Zane, these proceedings). 

Of course, if only part of the area that we see changed, 
then the average 'after' spectrum we currently observe con- 
tains a contribution from the unchanged parts of the surface, 
i.e., from the cooler, featureless 'before' spectrum . Hence, 
the spectrum from the changed part should be hotter and 
should have a stronger line than one would infer from the 
average. We can set an upper limit to the contribution from 
the 'before' spectrum by requiring that it does not exceed the 
'after' spectrum at any wavelength. From Fig. [6] one sees 
that the limit is about 70%, set by the 35^10 A region. 

The above could solve the second question: it might well 
be that after the change, the parts of the surface of J0720 that 
show an absorption feature in their spectrum, have a line as 
strong as that observed in J1308. It only appears weaker in 
the 'after' spectrum because it is diluted by the featureless 
emission from the unchanged parts of the surface. So, just 
a single question may be left: how can a neutron-star atmo- 
sphere, with presumably the same magnetic field strength 
and the same composition, and with only a modest, AT /T < 
0.2 temperature increase, emit such different spectra? 



5 Strongly magnetised atmospheres 

In interpreting the spectra, a major uncertainty is the com- 
position. For a single source, this may be difficult to deter- 
mine uniquely, but one can hope to make progress by treat- 
ing the INS as an ensemble: ideally, it should be possible 
to understand the features (or lack thereof) in all INS with 
a single composition, appealing only to differences in tem- 
perature and magnetic field strength (constrained by obser- 
vations where possible), which might lead to different ioni- 
sation states being dominant, and possibly the formation of 
molecules or even a condensate. Here, we discuss only the 
possibilities of hydrogen or helium atmospheres. For com- 
pleteness, we note that gaseous atmospheres composed of 
heavier elements appear to be excluded by the lack of large 
numbers of features. Condensated from heavier elements are 
also being considered seriously (Pons, Ho, these proceed- 
ings), and detailed theoretical calculations are being carried 
out to determine at what magnetic field strength condensates 
can form ( Medi n~& Lail20 06: Lai, these proceedings). 
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5.1 Hydrogen 

The presence of a hydrogen atmosphere has often been con- 
sidered by default, since if any hydrogen is present, grav- 
itational settling will ensure it floats to the surface. Typi- 
cally, it has been assumed the hydrogen is fully ionised, and 
spectral features have been interpreted as proton cyclotron 
lines. In strong magnetic fields, howe ver, the b inding en- 
ergies of atoms increase (for a review, Lai 2001; Potekhin, 
these proceedings; see also Fig.0, and for temperatures and 
fields appropriate for INS, a fr action of up to 10% of neu- 
tral hydrogen will be present (Potekh in et alJll9 99). From 
initial model-atmosphere calculati ons that ta ke the presence 
of neutral hydrogen into account ( iHo et al]l2003l see their 
Fig. 3), it is clear that, e.g., at 10 6 K and 10 13 G, the lines 
from neutral hydrogen have larger equivalent width than the 
proton cyclotron line (they are less deep but muc h wider, due 
to the so-called motional Stark effect; IPavlov & Meszarosl 
Il993t iPotekhin & Pavlovlfl997i) : at lower temperatures or 
stronger magnetic fields, the fraction of neutral hydrogen in- 
creases and hence the difference should be even larger. In 
general, it is worth stressing that the features are very strong: 
they may not appear so on the logarithmic scale typically 
used, but they have depths depth often exceeding 50%, sim- 
ilar to what is observed for J1308. 

Below, we first discuss whether the energies of the main 
features observed in the INS can be reproduced by a strongly 
magnetised atmosphere, and then turn to two possible prob- 
lems: harmonically spaced lines found recently, and the fea- 
tureless, black-body like spectra shown by some INS. We do 
not include the optical excess among these problems, since 
currently it is not clear any model makes reliable predictions 
for the optical emission: at a few 10 13 G, the plasma fre- 
quency exceeds the frequencies of optical photons, and the 
models do not take into account the resul ting significant de- 
viatio ns of the refractive index from unity Jvan Ad elsberg & I 
1200(1 see lKowalski & Saumonl2004l for a discussion of pos- 
sible effects in the context of cool white-dwarf atmospheres 
composed of helium). 

5.1.1 Line energies 

In Fig. we show the energies for features that might be 
produced in a hydrogen atmosphere: the electron and proton 
cyclotron lines, and the bound-bound and bound-free transi- 
tions of neutral hydrogen (relative to the ground state). Also 
shown are the approximate energies of the main features that 
have been detected in the various INS (corrected for an as- 
sumed gravitational redshift of 0.3) 

In the figure, thick vertical lines indicate the two mag- 
netic field strengths inferred from timing (§2). From those, it 
follows that if J0720 and J 1308 have hydrogen atmospheres, 
the features are most likely due to the transition from the 
ground state to the first excited tightly bound state of neu- 
tral hydrogen, perhaps in combination with the proton cy- 
clotron line. As argued in fj^Jabove, the line in J0720 might 
be weaker than that in J1308 because the emission from part 
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Magnetic Field (G) 



Fig. 7 Energy versus magnetic field for electron and proton cyclotron, 
ground state to tightly bound (tb) or weakly bound (wb) states in 
neutral hydrogen, hydrogen ionisation, and molecular H2 dissociation 
<Ho et alJI2051lLaill200HIPotekhinlfr^93) . For fields above ~ 10 I4 G 
(hatched region to the right), features may be washed out due to the 
effects of vacuum resonance mode conversion lLai & Ho 2003). The 
hatched bands show the energies of the main absorption features (cor- 
rected for gravitational redshift) for the INS, and the two vertical lines 
indicate the dipole field strengths determined by timing. The sources 
are labeled with short-hand notation a s in the text ( see end of §1). 
References for the e nergies are: J0420. fHaberl et all I2004al): J0720 
IFfaberletalJ i 2004b): J0806, Haberl et al]<2004ifr J1308. Hab erletalJ 
120031) : J1605. lvan Kerkwiik et al] 120041) : J2143. IZane"et alj 120051) 



of its surface is featureless (which in itself is problematic; we 
return to this below). Alternatively, the line in J0720 might 
be w eaker because it is to the second excited tightly bound 
-State (Ivan Kerkwijk et all2004l) . 

U If the above is correct, the feature in J0420 likely has the 
same origin and thus its field should also be a few 10 13 G. 
The features in J1605 and J0806 could result from the same 
transitions or from the ionisation edge, but in either case the 
implied m agnetic field strength is hi gher, close to 10 14 G. 
For J1605, van Kerkwiik et al]( l2004l) noted that the line was 
substantially weaker than that of J 1308, and they suggested 
this might be due to the effect of vacuum resonance mode 
conversion, which for fields in excess of -7 x 10 13 G tends 
to weaken features (see lHo & Lail2 003: also Lai, these pro- 
ceedings). Finally, for J2143, the line energy of 0.7 keV 
is substantially higher than what is observed for all other 
sources, and for any transition in neutral hydrogen, the up- 
per state is auto-ionising: it is at an energy level that is higher 
than the continuum from the ground state. It is thus not clear 
whether the line could be due to neutral hydrogen. Instead, 
it might be due to the proton cyclotron line in a field of just 
over 10 14 G. For these field strengths, the feature should 
be strongly weakened by vacuum resonanc e mode conver- 
sion (but not necessarily d isappear; e.g., IHo & LaH 120041: 

Ivan Adelsberg & Laill2006l) : qualitativel y, this is consisten t 
with the rather modest observed strength Z ane et all (12005) . 
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5.1.2 Possible problem 1: Harmonically spaced lines 

At the conference, evidence for harmonically spaced absorp- 
tion lines was presented for three INS. For J 1605, Haberl 
(these proceedi ngs) found that apart from the line at 0.40 keV 
discovered bv lvan Kerkwijk et alJ (l2004h . the EPIC-PN data 
show a significant feature at 0.78 keV, i.e., at an energy 
that is in a 1:2 ratio with that of the stronger line. Further- 
more, a third feature at 0.59 keV could be present, con- 
sistent with energies in a 2:3:4 ratio. For J1308, Schwope 
et al. (these proceedings) presented evidence that the sin- 
gle strong feature originally found at 0.3 keV or less by 
lHaberl et alJ d2003h . could be composed of two features, at 
0.23 and 0.46 keV, i.e., again harmonically spaced. Finally, 
for J0806, the single feature at 0.43 keV found by Hab erl et al. 
(2004a) may again be better described by two features at 
0.30 and 0.60 keV (Haberl, these proceedings). 

It would appear tempting to interpret these features as 
cyclotron lines, since those naturally have harmonic energy 
ratios. It is difficult, however, to see how this could be pos- 
sible for proton cyclotron lines, since the harmonics are ex- 
pected to be exceedingly weak: the oscillator strength for the 
harmonic would be a factor E /m p c 2 weaker than that for the 
fundamental. 

Instead, as mentioned in a discussion with George Pavlov, 
Joachim Trumper, and Frank Haberl at the meeting, a differ- 
ent solution may be suggested by the behaviour of the tran- 
sitions of neutral hydrogen. As can be seen in Fig. for 
any transition, above a certain magnetic field strength, the 
transition energy starts to become proportional to the pro- 
ton cyclotron energy. As a result, at sufficiently strong mag- 
netic field, the transitions become harmonically related. A 
possible problem, however, is that in this situation, the up- 
per level of the transition is an auto-ionising state, i.e., it 
has an energy in excess of the continuum energy relative to 
the ground state. It will still lead to some additional opac- 
ity, but at present it is not clear whether this is sufficient. 
Fortunately, there is one prediction: for J 1605, it would not 
be possible to explain the spectrum if there are really three 
features in a 2:3:4 ratio, without a strong corresponding T; 
thus, the prediction is that upon further analysis, the 0.59 keV 
feature will disappear. 

Finally, we note in this context that it will be worth check- 
ing carefully that for J2143, the 0.7 keV feature obse rved is 
in fact not a 'harmonic' From the present fits bv lZane et alJ 
J2005h . a rather high Nn is inferred, and this could perhaps 
be an artefact of a strong absorption feature at ~0.3 keV? 
(From initial attempts, this appears unlikely; Cropper, 2006, 
pers. comm.) 



ferred from the bow-shock shaped Ha nebula around the 
source (van Kerkwiik & Kulkarni 2001 ; Kaplan et al. 2002b), 
is below 10 13 G, in which case all features may be below the 
observed band (Fig.Ql, but for J0720 this explanation is not 
possible. Furthermore, for a mostly ionised atmosphere, the 
spectrum is expected to have a hard tail, unlike the observed 
exponential, Wien-like shape, since the free-free opacity de- 
creases with increasing energy. 

There are several possible solutions. First, there could be 
a reason for the opacity to be much greyer than currently es- 
timated, so that the emission at all wavelengths originates 
from layers at similar depths and thus with similar tempera- 
ture. The extreme version of this, discussed in detail by Pons 
and Ho (these proceedings), is that the two sources have a 
condensed surface (with possibly only a thin hydrogen layer 
on top). A less extreme version might be that the atmosphere 
does not contain just ionised and neutral hydrogen, but also 
molecules, which might have so many transitions that the 
opacity becomes effectively grey. Hydrogen molecules do 
indeed have a higher binding energy than hydrogen atoms, 
but the dissociation energy is only around 0.2 keV for a few 
10 13 G ( lLal200lt see Fig. [7}. With temperatures only a fac- 
tor two smaller, the abun dance should be very small (as in- 
deed found by, e.g., Potekhi n'et al.lll999l : see their Fig. 7). 
Nevertheless, it may be worthwhile verifying this, making 
sure that the abundance and the resulting opacity are indeed 
negligible. 

A possible alternative way to produce spectra resembling 
black bodies is by making the temperature profile in the at- 
mosphere shallower, closer to isothermal. While this cer- 
tainly appears ad hoc and likely would require significant 
fine-tuning, there is evidence for active magnetospheres: for 
J1856, the Ha nebula provides evidence of a pulsar wind, 
and for J0720, the optical emission appears to be partly non- 
thermal. If there is an active magnetosphere, some particles 
might hit the atmosphere, leading to additional heating; at 
the right location s, this could lead to r ather different emer- 
gent spectra (e.g.. lGansicke et al.l2002h . 

At present, none of the above explanations seem satisfac- 
tory. Also, none provide an easy explanation for why some 
sources have featureless spectra while others have not (or 
why it would change). Perhaps the first parameter to con- 
sider would be the overall temperature, since J 1856 is cool 
and the appearance of the absorption feature in J0720 was 
accompanied by a temperature increase. The increase was 
only small (AT /T < 0.2), however, and furthermore, an ab- 
sorption feature does app ear to be presen t in the coolest INS, 
J0420 (kT ~ 45 eV, Ha berl et all2004al) . 



5.1.3 Possible problem 2: Featureless black-body spectra 5.2 Helium 



Perhaps the most severe problem with the idea that the INS 
have pure hydrogen atmospheres is that the spectra of J 1856 
and J0720 (before the change) are featureless and well repre- 
sented by black-body emission. For J 1856, perhaps no fea- 
tures are expected, since its magnetic field strength, as in- 



Above, we stated that if any hydrogen were present, it would 
float to the top. Recently, it has been questione d, however, 
whether an outer hydrog en envelope can survive dChang & Bild sten 
2004; Chang et al. 2004). The reason this is not certain is 
that some hydrogen will diffuse down and reach underlying 
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Carbon or Oxygen layers , where, if the temp eratur e is right, 
it will be burned. Indeed, Chang & BildstenN2004h find that 
all of the hydrogen can be burned in the first 10 5 yr of a 
neutron star's life, in which case an atmosphere composed 
of helium might be left (unless hydrogen is replenished, as 
could happen due to spallation by relativistic particles from 
the magnetosphere, or very low levels of accretion). 

Pa rtly inspired by this possibilitv. lPavlov & Bezchastnovl 
d2005l> calculated properties of singly-ionised helium in strong 
magnetic fields. For a few 10 13 G, the transition energies are 
again in the range that features are observed in the INS, and 
hence it seems worthwhile to try to do a similar analysis as 
done above for hydrogen. From a very rough first attempt 
at producing model atmospheres (done by Kaya Mori and 
Wynn Ho), including neutral helium, it seems that, like for 
hydrogen, the features will be very strong. Typically, how- 
ever, more than one very strong feature should be present, 
which appears to be in conflict with what is observed. The 
picture is currently incomplete, however, since molecules 
have not yet been considered, while for helium the binding 
energy of, e.g., HeJ is sufficiently high that it may well be 
present (a detailed calculation is tricky, since one has to have 
a decent estimate of the number of possible rotational and 
vibrational states). 



6 Discussion and future prospects 

Of the four main parameters mentioned in the introduction 
that determine the properties of the thermal emission from 
INS, we now appear to have reasonable handles on three: 
the shapes of the X-ray spectra indicate temperatures around 
10 6 K, period derivatives imply magnetic field strengths of 
a few 10 13 G, and parallax measurements show that a fair 
fraction of the surface is emitting X-ray radiation. 

The main unknown appears to be the composition. We 
found that the energies of the observed absorption features 
can be matched fairly easily for hydrogen atmospheres. How- 
ever, reproducing the smooth, featureless spectra of some 
INS, and the Wien-like high-energy side of the X-ray spec- 
tra in general, appears problematic, nor is it clear how the 
spectrum of J0720 could change from featureless to one that 
has an absorption line. 

Fortunately, it should soon become clear whether these 
issues are real problems or not, since great progress is be- 
ing made in constructing more reliable strongly magnetised 
hydrogen model atmospheres (Lai, Potekhin, these proceed- 
ings). From ||5] it seems particularly important to include in 
full detail transitions to the auto-ionising levels, verify that 
all sources of opacity, including from (traces of) molecules 
are included, and check the influence, in particular on the 
temperature profile, of high-density effects and vacuum res- 
onance mode conversion. At the same time, it would seem 
worthwhile to consider atmospheres of other elements; for 
the INS, He might be most relevant, but it would be good 
to check heavier elements as well, since these may cause the 



absorption features seen in IE 1207.4—5209 (Hailev & Mori 
2002). 

From the observational side, the easiest route to further 
progress would appear to be timing. With further estimates 
of the magnetic fields, one can test the predictions based on 
hydrogen atmospheres, that J0420 has a field about as strong 
as that of J0720 and J 1308, J0806 a stronger one, approach- 
ing 10 14 G, and J2143 the strongest, in excess of 10 14 G. 

For the X-ray spectra, further monitoring is useful, but 
perhaps the largest advance will come from the unified anal- 
ysis of all sources, which allows one to exclude instrumental 
effects. This is already well underway for the EPIC-PN data 
(Haberl, these proceedings), and similar studies of the LETG 
and RGS data should prove fruitful. As present, first steps 
are being taken in detailed modelling of the phase-resolved 
spectra (Haberl, Zane, these proceedings), and this should 
help obtain stronger constraints on the thermal distribution 
over the surface. 

Finally, in the optical-ultraviolet regime, it would be good 
to complete the census of the sources, and obtain at least 
rough spectral energy distributions, to determine whether the 
emission is thermal, or whether there are non-thermal com- 
ponents. For sources that are sufficiently bright, proper mo- 
tion measurements can help determine true ages and parallax 
measurements can help determine distances. 
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